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By Eugene S . Love 

SUMMARY 


A reexamination has been made of the use of simple concepts for 
predicting the shape and location of detached shock waves . The results 
show that simple concepts and modifications of existing methods can yield 
good predictions for many nose shapes and for a wide range of Mach numbers . 


INTRODUCTION 


In recent years interest has arisen in the problem of predicting 
the form and location of detached shock waves . This interest has been 
stimulated by the necessity for blunt noses and leading edges on con- 
figurations designed for hypersonic flight in order to cope with aero- 
dynamic heating. The ability to predict the form and location of the 
detached shock is of primary importance in analyses of aerodynamic inter- 
ference and aerodynamic heating. Knowledge of form and location also has 
a more elementary use in that it often influences the choice of maximum 
model size for a given wind tunnel. Numerous experimental and theoreti- 
cal studies have been devoted to the determination of form and location 
as well as the important factors influencing form and location. (See 
refs. 1 to 36 , for example.) Many of the studies have been centered 
upon particular important details of the problem and, therefore, have 
been logically restricted in scope, for example, studies restricted to 
hypersonic speeds or to regions in close proximity to the nose . Other 
studies have been more general in that they present methods for calcu- 
lating detachment distance and shock shape without restrictions on speed 
or distance from the body. (See ref. 26, for example.) However, these 
methods us ually involve laborious characteristic calculations or such 
lengthy iterative procedures that they have been evaluated only for one 
or two conditions, and, although the results obtained are good, less 
laborious methods are preferable if the results obtained by the simpler 
methods are satisfactory. A number of methods have, in the light of 
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subsequent exp.erlmsn.tal results* been shown to be_ inadequate (see the 
evaluations made in refs. 24- and 25, for example) or severely restricted 
in application (see ref. 18 , for example). 


With the aid of the theoretical and experimental information now 
available, a reexamination "has been made of .the more successful simple, 
methods and concepts with a view toward (l) extending the range of appli- 
cability of— existing methods through modification* ( 2 ) presenting simple 
methods for .predicting shock shape' and detachment “distance, and ( 3 ) direct 
ing attention to areas where further study is needed. The results are 
presented in this paper. In the development of methods for prediction 
of— shock shape and detachment distance, a primary objective has been to 
obtain results sufficiently accurate for most practical engineering appli- 
cations . The methods presented are for air only (ratio of specific heats 
of 1.4) and are intended for use at supersonic and hypersonic speeds; 
the applicability of these methods at _ transonic speeds has not been exam- 
ined to any extent. 


SYMBOLS 


b distance between most forward point on shock and nose 

of body 

C constant defining detachment distance (see eq. (2)) 

C c value of _CL for convex face of complete hemisphere or 

two-dimensional semicircle symmetrically disposed with 
respect to free-stream direction 

Cq value of C for 8 ^ = S^ e t- C Q = 1) 

C^q value of— C for 5q = 90° (flat face) 

D diameter of sphere or cylinder 

d 1 diameter or height in plane of point on body determined 

by angle for shock detachment 

F horizontal distance measured from center of- spherical 

. (or circular) nose- or from face of flat nose to point 

on shock 

M ro . free-stream Mach number . 

q slope defined by equation ( 3 ) 
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x 
X ' 


y 


horizontal ordinate (free -stream direction) 

horizontal distance from most forward point on shock to 
plane containing d' 

distance from most forward point of detached shock to 
intercept of its asymptote on X-axis (see fig. 9)> 

vertical ordinate (normal to free -stream direction) 


p = \K 2 - 1 



ratio of specific heats 

semiapex angle of cone (or wedge) or equivalent cone 
(or wedge) 

semiangle of cone (or wedge) for shock detachment 

flow deflection for sonic flow immediately behind shock 

local inclination of detached shock measured with respect 
to X-axis that gives sonic velocity behind shock 


4 

9 


angle between normal to free -stream direction and control 
line (see fig. 9) 

, s. 

angle between X-axiS ahd line joining center of sphere 
(or cylinder) to corner of cut-off sphere (or cylinder) 

free -stream Mach angle 


* 


term defined by equations- (C7) and (C8) 


DISCUSSION 


Shock-detachment distance and its prediction will be cqnsidered 
prior to shock shape since the shape predictions hinge upon the detach- 
ment distance . 
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Shock -Detachment Distance 

Review of general concepts . - It is instructive to review at- the out- 
set several well-known and fundamental features pertinent to shock detach- 
ment and detachment distance. To this end, consider the simple cone- 
cylinder (or wedge-slab) having a sharp shoulder at the juncture of the 
cone and cylinder and with a semiapex angle of S Q and an attached shock. 

As 6 q is increased, the shock progresses continuously from a state of 

attachment to one of detachment in the manner described by Guderley 
(ref. ll) and by Busemann' (ref . 5)* As the value -of 8 q approaches 

that for detachment, a region of subsonic flow exists between the surface 
of the cone, the shock, and-the sonic line the origin of which muBt of 
necessity be at the shoulder (refs. 4 and 5 ) which is the center of expan- 
sion to supersonic flow. Clear ly, no p a nt of the bo dy downstream of the 
sharp shoul der affects the Setachment -process; a nd it becomes immediately 
apparent that the diameter at the shoulder' d 1 , and therefore- the diam- 
eter at the sonic point of the body, is one of the fundamental parameters 
in determining detachment distance . Dor the cone ^cylinder, the vertex of 
the detached shock for values of 8 Q barely larger than that producing 

detachment can be pushed no farther from the shoulder than the distance 
x 1 corresponding to the length of the cone that has a semiapex angle 
. Thus the maximum detachment distance (in the absence of viscous 

effects) is given by 

(r) - °-5 cot 6 det 

' 'max 

as has been previously inferred by Guderley (ref. 4) and others. Moeckel 
(ref. 9 ) refers to— this expression as the geometric method for predict ing 
detachment distance. 

Inasmuch as equation (l) gives the maximum detachment distance, 
x'/d' must be expected to decrease- as 8 q varies from to larger 

values. From a physical viewpoint it may be reasoned that, as 8 q 
increases beyond ^det' the tip of the cone recedes toward the shoulder 

and loses contact with the shock; once- contact- is lost? and as the cone 
tip recedes, the shape of— the shock in the immediate vicinity of the nose 
adjusts itself to a shape of more uniform curvature in such a way that 
the peak; in the detached shock that _ ~occurs when the cone tip is in close 
proximity to the detached shock (see- -shadowgraphs _of ref. 27 ) is gradu- 
ally eliminated and the front of the shock moves closer to the shoulder. 
The observations of Busemann (ref. 5 } ' follow readily: there is no appre- 

ciable effect of nose shape upon detachment- distance measured from the 
shoulder unless the nose tip is in close proximity to the detached shock; 
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and , except for the case of close proximity, the detachment, distance and 
the' shape of the detached shock between the sonic -points on the shock 
are determined by the shoulder of the body . _In order to locate the effec- 
tive^houlders of bodies that - have rounded noses, that is, bodies not 
having sharp shoulders of the type common to the cone -cylinder, Busemann 
proposes the use of the most upstream point on the body surface that is 
tangent to a line inclined at S^ ei _ with respect to free-stream direc- 
tion but also recognizes that, although this is the important point of 
the body, it is not in this case the location of the sonic point. This 
concept affords a simple means of correlating detachment distance for 
bodies of various shapes, and its general adequacy has been substantiated 
by experimental results (refs. 9 anti 11, for example) for several dif- 
ferent nose shapes and for Mach numbers up to about J. 

Compilation of and general correlation of data on detachment 
distance . - In the past few years, additional experimental information 
on detachment distance has been obtained and it is of interest to see 
whether these additional data may be correlated by the principle of 
Busemann. A compilation of experimental data according to this principle 
(that is, in terms of x'/d') is presented in figure 1. When these data 
were compiled, it was observed in several sources that some of the experi- 
mental data gave values of x'/d' that fell above the curve given by 
equation (l) (for example, some of the data of refs. 28 and 35)- Since 
these points are believed to represent improbable values and were in a few 
instances recognized by the investigators as being associated with large 
experimental inaccuracies (ref. 33* for example), they have been omitted 

from figure 1. Where necessary, values of 5 for converting the 

Q.6t 

detachment distance to the form of x'/d' were obtained from reference 37- 

The data of figure 1 show that the parameter x'/d' is a unifying 
one and is justified throughout the Mach number range of the experimental 
data shown. Similar conclusions of earlier studies in the supersonic 
speed range thus remain supported into the hypersonic speed range. 

Although the experimental data tend toward a single band of data for both 
the two-dimensional and axisymmetric compilations, there are differences 
at a given Mach number that can only be attributed to the effects of nose 
shape; examples of these "effects and their prediction are shown subse- 
quently. The continuity method of Moeckel (ref. 9) is seen to give a 
prediction that is in generally good agreement with the band of experi- 
mental data. The experimental data shown in figure 1 are for blunt nose 
shapes or for _ conditions where is appreciably larger than 6^^. 

For spheres and for circular cylinders normal to the stream (two- 
dimensional) , there is a geometric minimum x'/d' defined by the nose 
of the sphere and of the circular ■ cylinder above which all experimental 
data must obviously lie. The curves defining these minimum values are 
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shown at the bottom of figure 1 in order to convey some idea of-the prox- 
imity of the shock to the nose of these shapes . 

Refinements for particular nose shapes . - On the basis of the avail- 
able experimental, data, certain refinements can be made that will enable 
a more accurate prediction to be made of detachment distance for certain 
nose shapes. A convenient basis for refinement is to modify the expres- 
sion for the maximum nondimens ional detachment distance as given by 
equation (l) to the following expression for the detachment distance 
under any condition: 


fr - °- 5 ° cot 8 det (2) 

In this expression the factor C is to be determined. The upper limit 
of C is obviously 1 ^this limit will be designated C Q ^, and the 

lower l imi t will not be greatly removed from 1 in view of the secondary 
effects of .nose shape and of the proximity to the curves for maximum 
detachment distance shown by the experimental data in figure 1. 

Value of C for flat faces .- The value of C for flat faces will 
be designated as C^q. For axisymmetric bodies with flat - faces, such 

as circular disks normal ter the stream or the cone -cylinder with 

= 90°, a value of C^q = 0.70 appears to give good agreement with 

experimental results over the Mach number range of the experimental data 
contained in figure 1 . For infinite Mach number this simple conversion 
yields x'/d' = 0.222; this value compares favorably with a prediction 
of Serbia (ref-. 16) for a circular disk normal to the stream and, when 
expressed in terms of the present parameters, yields 0.2J0 (shown to the 
far right in fig. 1 ) . There is not as much experimental information on 
two-dimensional flat faces as for axisymmetric ones but, for the avail- 
able two-dimensional data, a value of - C^q = 0.86 appears to be satis- 
factory. (See fig. 1 .) For an infinite Mach number this value yields 
x'/d’ = 0.421 for a two-dimensional flat face. 

Value of C for circular faces . - The value of C for the convex 
face the generatrix of which is the l80° (or greater) arc of a circle 
symmetrically disposed with respect to free-stream direction will be 
designated by C c (the subscript c meaning complete to distinguish 

complete circular faces from modified circular faces to be considered 
subsequently). For complete circular faces, both axisymmetric and two 
dimensional, the value of C appears to vary significantly with Mach 
number. Figure . 2 presents the variation of C c with Mach number for 

the axisymmetric circular face (or sphere) and for the two-dimensional 
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circular face (or the circular cylinder normal to the stream) . These 
curves were determined from the appropriate experimental data contained 
in figure 1. It will be noted that the curves have been extrapolated 
slightly beyond the range of the data contained in figure 1. 


An interesting comparison can be made between the value of C c for 

a sphere at hypersonic speeds as indicated by the curve of figure 2 and 
the values for detac hm ent distance which Serbin (ref. 3 6 ) and Hayes 
(ref. 21) have calculated for a sphere for the case of approaching 

°°. The results of Serbin and Hayes may be expressed in terms of C c . 

When this is done (and with y = 1.4), Hayes' prediction yields 

C c = 0.825 and. Serbia's prediction gives C c = 0 . 850 ^ the latter value 

is in close agreement with the value C c = 0.857 which the curve of 
figure 2 approaches at hypersonic speeds. 


Variation for cone -cylinders and wedge-slabs .- For cone -cylinders 
having a semiapex angle 5 q between S^ e ^_ and 90 ° > the refinement is 


not as readily obtained, but the recognition of certain features permits 
a prediction to be made of the effect of varying 5 q between 8^ ^ and 

90°, and Johnston's experimental measurements of this effect (ref. 27) 
afford a ready means of checking a prediction. Consider first the vari- 


ation of 



replaced by S Q 


with 

x' 

and — 
d' 


5 det as given by ^equation (l) . 
is replaced by the cone length 


If 

x 

D* 


6 det is 
an expression 


for the variation in cone length with 8 q is obtained. This variation 
is represented by the curve AA' in the following sketch: 



If the semiapex angle 8 q of a cone -cylinder in a flow at constant 
supersonic Mach number is allowed to increase continuously from some 
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very small value to 90 ° , the following conditions occur. At first the 


shock Is attached and the distance — between the shoulders of-the cone- 

D 


cylinder and the vertex of the attached shock varies according to the 
curve AA‘. However, when Sq reaches S^ e j__ (point D) corresponding 

to the given Mach, number, the distance between the vertex of the- shock 
and the shoulder ceases to vary according to the curve AA'. The value 


of 




for a detached shock (point E) has been reached and, with 


max 


a further increase in Sq, — j- must decrease according to some curve. 


say BB', which defines the variation in detachment distance. It is 
apparent that an infinite number of curves of the type BB' may branch 
off from the curve AA', since for every value of M OT there is a corre- 


sponding curve BB'. Point B will never reach A.V since for M m = » 
the value of 8^^ is 57*6° for cones (after Maccoll, ref. 35) and about 
for wedges . The- problem is thus one of determin,ing the general 


45.6° 


form of the curve BB* . 


Guderley (ref. 4) has shown that the transition from attached shock 
to detached shock is one of continuous change. Therefore, it is reason- 
able to assume that at point B the slope of the curve BB' will be 
equal to the slope of the curve AA'. This slope will be designated by 
q and, from equation (1), it has the value 


q = -0.5 csc 2 5 det (3) 

Thus the ordinates and slope of the curve BB ! at point B can be cal- 
culated. When 6^ = 90° (point B'), the detachment distance x'/d* 

corresponds to the flat-face condition for which empirical values of C 
in equation (2) have already been proposed. The slope of the curve BB' 
at B' must be essentially zero, if not-exactiy so. Thus the ordinates 
and slope at both ends of the curve BB’ are obtainable. Prom the works 
of Guderley and Busemann, the form of the curve BB' must be such that 
x’/d' is always decreasing and the rate of decrease in x'/d' with Sq 

is greatest at ®det‘ Consequently, from the knowledge- of the end-point 

conditions and the restrictions on the form of the curve, one is led to 
suggest that the curye BB' is close to being elliptic in form. The 
development of the" general elliptic equation giving the variation in 
x'/d' between 5 q =• S de -fc and 8 q = 90° is given in appendix A. 

Predictions of the "variation in detachment distance for cone- 
cylinders and wedge-slabs (two-dimensional) according to, appendix A 
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are given in figure 3 for several Mach numbers , and comparisons are made 
between Johnston’s experimental results (ref. 27) at = 2.45 and the 
predicted variations for this Mach number . The predicted”" variations 
agree closely with the experimental results. The two-dimensional wedge 
data fall slightly above the predicted curve; however, one suspects that 
this difference may be attributed to the difficulty in obtaining two- 
dimensionality in experimental investigations of this type (see refs. 11, 
24, 25, and 2 6, for example) since the experimental data for the wedge 


do not show the correct values of 


m 


for 8. 


max 


det 


at M„ = 2.45* Fur- 


ther, the experimental data for the wedge are, in general, uniformly higher 
than the predicted curve by the amount of the experimental error at 


The same procedure as employed for the prediction of figure 3 may 
be used to calculate the variation in x'/d' with M for constant 6 ,.. 

Examples of this prediction are given in figure 4 for a cone -cylinder 
at several values of 8^ . No suitable experimental data for cone- 

cylinders were found for comparison; however, the results of figure 3 
lend validity to the predictions of figure 4, and the experimental 
results of Griffith (ref. 28) for wedges tend to substantiate the type 
of variation shown here when M ra is decreased below that for detachment. 

Variation for cut spheres and cylinders .- One interesting model 
that has been "used in the study of detachment distance is the so-called 
cut sphere (axisymmetric case) or cut cylinder (two-dimensional case). 

The cut sphere or cylinder is so termed because the shape of the sphere 
or cylinder is altered by actually cutting segments from the sphere or 
cylinder or the models are so constructed that they simulate the effect 
of this cutting. For example, the diameter of a hemisphere -cylinder 
may be reduced by concentric, machining such that the nose shape varies 
systematically from a hemisphere to a flat face (of zero diameter in 
this limit, however) . The nondimens ional results thus obtained (with 
the exclusion of the zero -diameter limit) can be considered to be the 
same as those that would be obtained, for example, with a number of 
models of constant diameter and with varying radius of the nose. A 
similar procedure is applicable to the circular cylinder normal to the 
stream (two-dimensional) which in the case of actual cutting amounts to 
removing symmetrically disposed segments along parallel planes that are 
also parallel to the plane of symmetry and to the free -stream direction. 

The cut sphere or cylinder' is thus seen to afford a convenient 
means for examining the effect upon detachment distance of systematically 
varying the radius of the nose while diameter or height is held constant, 
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and thus a range of Interesting nose shapes is covered. Further, one is 
able to determine the radial point on a sphere or cylinder that is most 
important in fixing detachment distance. 

In reference 24, Kim reports the results of two -dimensional- tests 
of a cut cylinder at = b i However, in reference 2b there are some 

discrepancies in the figures and in the values tabulated and those indi- 
cated in the figures of that reference. In order to justify the correc- 
tions that are made herein to account for the discrepancies appearing 
in reference 2b, the shock shapes as reported by Kim are reproduced 
herein in figure 5- In order to agree with Kim’s tabulated values, the 
abscissa scale would need to be condensed as shown. With this modifica- 
tion a check is obtained of Kim's quoted values exceplrfor 9 = 50°, 
for which a value of 0.20 is indicated. (See tabulation at top of 
figure 5») Kim's results are expressed in the-ratio of the distance 
between the shock and the nose to - the diameter of the complete or uncut 
cylinder (9 = 90°) . Although this ratio is an irrelevant one from an 
analytic viewpoint, its use is more appropriate here- than the use of the 
correlating form x'/d', since the object is to determine the critical 
value of- 9 which in essence is determining d'. (The critical value 
of 9 defines t he- po int on the nose that— s eparates t he portion of t^ e 
nose~~th a t. Rf’^ct. s .detachmerviPdls tance from the portion that does not - ;) 
^further, by use of Busemann's concepl of — d' and sTTew simple asBUn$- 
tions, it is possible to calculate the variation of b/D with 9. 

In order to calculate this variation within the framework of the previous 
methods employed herein, one recognizes at the outset that the value of 
C fa r t .he nmcut cylinder (otT sphere) (8 = 90°) is different rronT^thg" 
val ue of C for the flat - face (9 = 0 U )T T he variation m. u 'betwe en 
these two limits i.s thus needed. As a first approximation it may be 
'assumed that this variation is linear; from this assumption it~ follows 
that 


C c “ C 90 + 9 


(4) 


where C c is the value of- - C for the complete— cylinder (or sphere), 

dC = C c - C^q, C^q is the value of C for the flat face, and 

d9 = 90° - 5^ e ^_ . From the previously presented values of C c and C^q, 

the linear variation of C with 9 is easily obtained. With the aid 
of equation (2), the variation of x'/d' with 9 may be calculated and 
converted to terms of b/D by the following equation: 

£ = sin 9 + S ^!- 9 - i 
D d 1 2 2 


(5) 
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where 


0° ^ 0 < 90° - 6 


'let* 


For values of 9 between 90 and 

90 u - 8^ e .j., the present analysis gives constant values of b/D at the 
value calculated for 0 = 90° - 8 det - (Note that, for this statement 
and for equation ( 5 ), D is the diameter of the uncut model.) 


From the results of figure 3> an elliptic variation of C would 
appear to be a more justifiable and accurate representation than the 
linear one. The method for computing the elliptic variation is covered 
in appendix B. 

Figure 6 presents the calculated linear and elliptic variations of 
C for the two-dimensional cylinder of Kim's experiments at M ro = 4. 

Also shown for comparison is the variation in C for wedge-slabs 
between ®det and 8 q = 90° as considered previously. For con- 

venience, both 0 and Sq are shown as the abscissa, 9 being merely 
90° - 8 0 . 


Figure 7 presents Kim's experimental measurements and compares them 
with the present predictions. The correction to the experimental point 
at 0 = 30 ° as shown by figure 5 i8 included in figure 7> as is a cor- 
rection to the point at 9 = 42° which is apparently misplotted in a 
similar figure in reference 24 since it does not agree with the quoted 
value in reference 24 or the value obtained from the shock locations / 
(reproduced in fig. 5 herein) . Of the present predictions the linear 
variation of C gives the better agreement with Kim's data fairing but 
this seems to be meaningless for the reasons just stated. When the data 
are plotted according to what are believed to be the correct values, the 
elliptic prediction is to be preferred and probably lies within the 
accuracy of the experimental data. Kim justifies the critical value of 
0 = 48°, obtained by extrapolation, on the grounds that the critical 
value of 0 must coincide with "the foot of the last Mach line" or 
sonic line, and in support of this statement quotes a calculated critical 
value of 0 of 44°. However, in view of the fact that the value of 
0 = 48° was , obtained by extrapolation of a faired curve through appar- 
ently erroneous points, as shown in figure 7 herein, this argument is 
weakened. The present predictions place the critical value of 0 at 
51. 2° ^which is 90° - 5^ e ^. as dictated by Busemann's concept). It is 

not only possible but highly probable that the foot of the sonic line 
occurs at 9 < 51*2° for, as Busemann has pointed out in reference 5, 
a peculiarity of the detachment phenomena Is • that the critical point 
which plays the primary role In fixing detachment distance for arbi- 
trarily rounded nose shapes turns out to be situated where the flow is 
purely supersonic. Both the foot of the sonic line and the foot of the 
last characteristic which intersects the sonic line are ahead of this 
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critical point. However, as also pointed out in reference "the body- 
slope mast-decrease to— smaller angles than _ if the sonic line is 

to complete its enclosure of the subsonic region "behind the shock. Thus 
what at first glance appears to be a peculiarity actually is a fulfillment 
of a necessary condition. 


One may reasonably inquire as to the difference that there would 
have been in the predictions of b/D in figure 7 if no account had been 
made for the change in C from C^q to C c . A brief examination shows 

that the predictions would suffer considerably. For example, if C had 
been held constant at b/D would have been 0.23 at 0 = 51*2° 

instead of 0 . 269 . 


Recently, Mr. Robert V. Rainey of the Langley Laboratory has 
obtained results (unpublished) for a series of axisymmetric shapes at 
Moo = 3.55 in which. the diameter was held constant— and the radius of 


the nose was varied in such a way that the nose shapes varied-from a 
hemisphere to a flat face . These results are particularly interesting 
since they include data for a value of 0 very close to the critical 
value predicted by the Busemann concept ^90° - . The data are again 

presented in terms of b/D rather than x'/d' in order to bring out an 
important feature of this type- of presentation. These results are shown 
in figure 8 and are compared with the predicted curve employing an ellip- 
tic variation of C determined according to appendix B. Values of 
x'/d' thus obtained may be converted to — b/D values by the relation 


i + 1 2 

D d' 2 tan 0 2 sin 0 

for 0° < 0 ^ 90° - 5, . and by the relation 

det 


( 6 ) 


b 

D 


Sin - S det) 

sin 0 


2LL + «L_ 

d’ 2 tan ( 90 ° 


®det) 


1 

2 sin 0 


(7) 


for 90° - ^ 0 = 90° • (Note that, for these equations and for the 

experimental models of fig. 8 , D is effectively the diameter of cut 
models as contrasted to the results of fig. 7; thus, b/D is not con- 
stant for values of 0 between 90 ° and 90 ° - as in fig. 7 *) 
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The prediction shows excellent agreement with the experimental 
results, and the experimental results confirm the existence of a pre- 
dicted knee at 0 = 90° - in this case at 38.5°. The Buse mann 

concept of the most important point of an arbitrary blunt profile bein g 
determined by ^det is once again clearly substantiated. The results 

of figure 8 should serve to correct the impression left by less complete 
experiments that the variation of b/D with 9 between 0° and 90° 
occurs smoothly and without a knee near 6 = 90 ° - 5^ . 


Shape of Detached Shocks 

Initial Considerations . - When the methods that have been proposed 
for predicting Shock shape which do not involve laborious procedures 
were examined, the method of Moeckel (ref. 9) appeared to offer the best 
possibility, when modified, for giving satisfactory predictions at both 
supersonic and hypersonic speeds for nose shapes that do not approach 
too closely the condition of attached shock. This method has been shown 
to give generally satisfactory results at supersonic speeds (see refs. 9 
and 11, for example) but Is known to be inadequate in Its present form 
at hypersonic speeds . The method is proposed only for that region of 
the shock between the sonic points; however, Moeckel 1 s experimental 
results show that the method may for many nose shapes be satisfactory 
as an approximation of the shock shape at distances considerably beyond 
the sonic points on the shock. 

Before modifications are discussed, it would perhaps be worthwhile 
to review briefly Moeckel' s method. (See ref. 9.) In this method it 
is assumed that’ the shock shape is hyperbolic in form and that its loca- 
tion and scale in relation to the body are determined from continuity 
considerations involving certain assumptions . These assumptions include 
the location of the sonic point on the body by the use of 5^^, the use 

of a straight sonic line, and the determination of the sonic point on 
the shock by assuming the straight sonic line to be inclined at an 
angle 7 with respect to the vertical where q is equal to 

|>(®det + 5 s) 6 s i s the deflection of the flow at the sonic 

point on the shock. As has been recognized by Moeckel, two severe assump- 
tions of his method are the use of a straight sonic line and the sonic 
line being inclined at l/ 2 ^ 8 ^ e ^_ + 8 g ^ . 

Development of present method . - In the present development, major 
modifications are made to Moeckel ' s method. The continuity relation is 
not used, -and only the trigonometric derivations are retained. The pres- 
ent method employs the results of the first part of this paper to cal- 
culate detachment distance instead of the continuity relation of 
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reference 9* Secondly, the straight line— arising from the point - on 
the body determined by 5^^ (see fig. 9) and- determining the point on 

the shock where the shock inclination e is equal to- that- for sonic 

s 

velocity behind the shock is not restricted in inclination to a value of 
T) = 1/2^8^^ + 8^ . Further, this straight line is not regarded as being 

a sonic line but as a control line that reproduces the effects of the 
true sonic line. There are several reasons for this assumption. The 
sonic point on the body lies ahead of that given by (except- for 

the case of the sharp shoulder) , and the exact sonic line is curved and 
according to reference -26 is, as a general rule, not - normal to the stream 
lines at- the sonic points on the body and at the shock; neither is it 
necessarily normal to the streamlines between these-points . It is clear 
then that a straight control line- which reproduces the effects of the 
true sonic line will not satisfy the inclination of the true sonic line 
at both of its extremities nor, as a general rule, is the control line 
likely to be inclined at the mean of the inclinations at the two extrem- 
ities of the true sonic line. Further, the inclination T] of the con- 
trol line that most effectively reproduces the influence of the true 
sonic line will not necessarily be that which gives a line that appears 
to represent best the actual location or average inclination of the true 
sonic line . It follows that - the determination of the value of- t] for 
the control line by analytic methods would be difficult - ; However, 
since the publication of Moeckel's work (ref. 9), a large number 
of experimental shock shapes have been obtained over a wide range of 
Mach number. With these experimental results and the advantage of hind- 
sight, together with an expression for q derived from Moeckel's trig- 
onometric relations , one is able to determine values of r| from known 
shock shapes. This expression for— rj is given in appendix C, as are 
simple relations for converting shock ordinates for circular and flat 
noses. Several calculations of T} quickly revealed that- the hyperbolic 
form of shock shape is best adapted to the spherical or circ ular nose 
(rather than the flat) when it is desired to obtain a value of f] that 
is suitable for shock shape both near to and far from the nose. Conse- 
quently, attention was centered upon obtaining the variation in q for 
spherical and circular (two-dimensional) noses. In order to determine 
a value of tj for a given Mach number, calculations were made for several 
points along the shock and the average was taken of the resulting values of 
T). In general, the calculated values of r| ranged approximately ±6° 
about the average value at the higher Mach numbers with less scatter air 
the lower Mach numbers. It is important to note that, whereas the value 
of T] is obtained by means of the tangent function (see eq. (C3) of 
appendix. C) , and therefore varies within the limits of- the principal 
values of the tangent (± 90 °), the actual movement of the control line 
(fig. 9J Is such. that rotates continuously counterclockwise with 
increasing M^. Thus negative angles .(calculated from eq. (C3)) have 


J 
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been converted to their positive complements (> 90 ° but < 180°) to achieve 
compatibility with the rotation of the control line. 


Variation of t) with M^.- The variation of obtained by the 
above procedure is shown in figure 10 for the sphere and the circular 
cylinder (two-dimensional) together with the values of ®det* S s' and 

Moeckel's value of l/ 2 ( 5 £ e t + 8 s ) . The adequacy of the values of r\ 

obtained in the present analysis will be shown subsequently by using 
these values in the calculations of shock shapes and comparing these 
calculated shapes with a number of the experimental shapes. The experi- 
mental shapes will include some of those from which points were taken 
to calculate the average n -values . The striking feature of figure 10 
is the increasingly large difference between Moeckel's values (given by 
T) = l/2(8 det + S s )) and those of the present analysis. For spheres, 
the present analysis gives values of tj that are always greater than 
l/2(8 det + 8 S ); near Moo = 8 the present values are about 2^ times 

greater. For two-dimensional circular cylinders, the present values 
fall below l/2^5 de ^ + & s ) at the lower Mach numbers and above, at the 

higher Mach numbers. With regard to these comparisons with Moeckel's 
results, it should be recalled that his method is proposed only for the 
region between the sonic points of the shock, whereas the present analy- 
sis attempts to include the shock beyond the sonic points as well. There 
was some indication that the values of t} determined in the present 
analysis would be smaller than those shown if the analysis had been con- 
fined to the portion between the sonic points; however, these smaller 
values would still be much larger than l/2 ( 5 de -^ + 8 S ) at the higher 
Mach numbers . 


Shown in the upper right-hand part of figures 10(a) and 10(b) are 
three values of T) that may be indicative of the magnitude to be 
expected at infinite Mach number. The value designated tangential r\ 
corresponds to the condition for which the control line becomes parallel 
to the nose surface at its point of origin on the surface (that is, 

90° + 8 de ^.) . The value designated as being determined from C c corre- 
sponds to the condition for which the control line intersects the axis 
of symmetry at the vertex of the shock whose detachment distance is 
determined from the value that C c tends to obtain at hypersonic speeds 

in figure 2 (c c = 0.857 for spheres and_Jl,_ 9_52 for cylinder s) ; this con- 
dition implies that the sonic point on the shock is at the vertex of the 
shock. The value designated maximum T) corresponds to the condition of 
the control line passing through the tip of the nose on the axis of 
symmetry; this condition implies that the shock touches the nose, that 
the sonic point on the shock is at the vertex of the shock, and therefore 
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that the sonic points on the shock and on the nose coincide. The con- 
ditions associated with maximum q and t] determined from C c require 

the control line to pass through the body; this passage through the body 
is not a point for concern since the control line was proposed as repro- 
ducing the effect of the sonic line and not as simulating its location 
or its termini. 


Some interesting similarities exist between the variation in t), 
or more properly, the control-line termini obtained in the present anal- 
ysis and the termini of the sonic line- obtained by Serbin (ref. 36 ), 
Chester (ref. 31) , and apparently by Freeman (ref. 32). Chester's anal- 
ysis shows explicitly that for all values of 7 the position of the 
sonic point on the body does not- experience large changes as 

increases from moderate supersonic speeds to speeds approaching infinity, 
whereas for all values of 7 except unity the sonic point on the shock 
experiences significant changes and moves toward the vertex of the shock. 
(For 7 = 1 , the sonic point-on the shock is indicated to be at the vertex 
of the shock.) Thus, for 7 other than unity, these variations of- the 
sonic points are of the same type as the variations exhibited by the 
termini of the control-line in the present— analysis . Both Serbin' s and 
Chester's analyses indicate that, when = °o and 7 = 1 , the sonic 
point on the shock is at the vertex of the shock, the detachment distance 
is zero, and the sonic points on the shock and on the body coincide (same 
sonic point conditions as discussed in the preceding paragraph for maxi- 
mum ti). For = co and 7 = 1.4, the same analyses indicate a finite 
detachment distance and the sonic point on the shock and on the body to 
be removed from the axis of symmetry. For comparison with the control 
line inclination of the present analysis, a straight line connecting the 
sonic points of Chester's analysis would be inclined at 90° with respect 
to the vertical. Thus, although some interesting similarities are 
observed between the locations and movements of the sonic points of these 
analyses and the behavior of the control-line inclination and termini 
of the present analysis, there are also differences that remind one not 
to lose sight of the major differences between the true-sonic line and 
the control line, as pointed out-previous ly. 


Effect of C and q upon shock shape .- Ihe procedure for calcu- 
lating the shock shape is outlined in appendix D. It is of interest- to _ 
examine at the outset the effects of changes in C and- T) upon shock 
shape. This is done in figures 11 and 12 by comparing the various pre- 
dictions with experimental results at M ro = 1.94 for a hemisphere - 


cylinder. Figure 11 presents the shock shape in terms of d' and 
with horizontal distance referenced to the vertex of the shock in the 

~ (See fig. 90 

cl a. 


This is the same form employed by Moeckel 


form 
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in his evaluation of his predictions in references 9 and 11. Figure 12 
presents the shock shape in terms of the maximum diameter D and with the 
horizontal distance referenced to the center of the sphere F/D. 

The results of figure 11 show readily (as do the results of fig. 12 
by closer examination) that, although C was developed as a means for 
obtaining detachment distance, it also has a significant influence on 
shock shape. Figure 11 also shows that, when this form of presentation 
is used, one may obtain an almost identical prediction with different 
combinations of T) and C and that these combinations would give an 
excellent prediction. of the experimental shock shape (not to be confused 
with shock location) . One may immediately conclude that, although this 
form of presentation is convenient for isolating- the effects of C, it 
could lead. to ambiguities and thus is not in itself always adequate in 
the present analysis to evaluate the predictions of shock shape. Accord- 
ingly, the form of presentation shown in figure 12 has been used. In 
this figure are shown the effects of C upon shock shape and location 
for q = l/2(6 det + 8 fl ) within the C-liml.ts established herein (l.O 

and 0-7) for axisymmetric noses and the effects of q upon the shock 
shape for the value of C c = 0 . 804 determined herein for spheres, at 

M m = 1.94. (See fig. 2.) With C c = 0.804, the position of the shock 

at the nose is, as is to be expected, accurately given. With increasing 
q (from 5 s , to Moeckel's value, to S^Let* arL(i finally to the value 

given by fig. 10), the shock curves downstream more rapidly and the pre- 
diction given by the present analysis (t) = 44° from fig. 10) Is in 
excellent agreement with the experimental results . The prediction given 
by q = is also good at this Mach number and the prediction given 

by q = l/ 2 (S<iet + S s) ^- s f alr - From these results and from figure 10, 
one sees some basis In the past practice of replacing 
by when Moeckel's method is used; figure 10 shows that &^ e ^. 

and the curve of the present analysis are in close agreement below about 
M ot = 2, and above this Mach number the curve for is always nearer 

the curve of the present analysis. However, at high Mach numbers this 
is indeed a trivial point. One also sees why the use of Moeckel's value 
gives reasonable predictions at low Mach numbers, particularly over that 
portion of the shock near the nose for which it was intended. 

Examples of shock prediction by present method .- Numerous calcula- 
tions have been made by the present method of the location and shape of 
detached shocks . These calculations have been compared with experimental 
results, and some additional random evaluations of changing q have been 
made. Some examples of these comparisons and evaluations are shown in 
figures 15 to 17 - 
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Figure 13 reproduces the shocks obtained by Kim (ref. 24) on a 
two-dimensional circular cylinder at - Mach numbers from 1-35 "to 6 . Shown 
for comparison are the predicted shocks with t} = §&et' which in essence 
equals Moeckel's value of l/2^5 det + 6 g ) for two-dimensional flow, and 

with r\ given by figure 10(b). With T) = 5 det or l/2(6 det + 8 g ) , the 
predicted shocks curve rearward too rapidly at low M*, and too slowly at 
high Mo,. For example, examine shock shapes at- = 1.8 and = 6 . 
With Tj varying according to figure 10 (b), the predicted shocks are in 
much closer agreement -with the experimental shocks at all Mach numbers . 
Figure 15 (a) also contains two-dimensional results and shows excellent 
agreement between prediction and experiment. 

Figures 14(b), 15 (b), 16(a), and 17 compare the predicted shocks 
with the experimental shocks for hemisphere cylinders for Mach numbers 
of 3*55> 5*8, 6 . 8 , and 7* 7* respectively. In all instances the present 
prediction agrees well with experiment. Figure 16 (a) also includes the 
prediction with x\ = 8 de ^_ to demonstrate the increasingly large dis- 
agreement with experiment that accompanies the use of this value of T} 
with increasing a prediction with q = l /2 (5 dei - + would show 

larger disagreement with experiment. 

Figures 14(a) and 16 (b) show what“may be expected in the way of 
predicting the shock for flat -face cylinders alined with the flow at 
Ma, = 3-55 and 6 . 8 , respectively, by converting the corresponding pre- 
dicted shocks for a hemispherical nose (figs. 14(b) and 16 (a)) by use 
of the Busemann concept. For the flat -face condition d' becomes D 
and coincides with the location of the flat face. Thus, the necessary 
conditions for conversion are simply 


and 


( 5 ). 


flat 



spherical or circular 


( 8 a) 


/F\ = /jl - - *1) 

' D 'flat ^ d 'spherical or circular 


( 8 b) 


The resulting shock shapes for the flat face obtained by this 
conversion are not in good agreement - with the experimental shocks, 
particularly at-large distances from the nose. However, some disagree- 
ment is to be expected. The values of t} for a flat face would prob- 
ably be less than those for a round nose in view of the influence of 
the sharp shoulder of the flat face in defining the area between the 
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shock, and the nose where the flow is choked. Further, the disagreement 
at large distances from the nose is in the direction to he expected since 
the sharp shoulder of the flat face is the source of a centered expansion 
which interacts with the shock. in a manner that reduces its inclination 
more rapidly than would the gradual expansion from the spherical (or cir- 
cular) nose. There is also the possibility that the sharp shoulder pro- 
duces a separation bubble and an associated overexpansion, but experi- 
mental evidence on this point is at present not sufficient to indicate 
its significance at moderate or high Mach numbers. In spite of these 
shortcomings, it would appear reasonable to regard this prediction as 
a suitable first-order prediction for flat faces. The important point 
to be gathered from these flat-face predictions is that they indicate 
the maximum discrepancy that may be expected in predicting shock shape 
by means of the spherical (or circular) calculation for noses that vary 
all the way from hemispherical to flat shapes (in the manner of the nose 
shapes of fig. 8) . As the nose shape departs from the flat and moves 
toward the hemispherical, the agreement between prediction and experi- 
ment becomes increasingly good. 

Figure 17 compares the present prediction with the experimental 
results at = 7*7 presented by Lees and Kubota in reference 22 for 
a hemisphere -cylinder . Also shown are the predictions made in ref- 
erence 22 by means of blast-wave theory. The present prediction is in 
excellent agreement with the experimental shock. 

Shock shapes for blunted cones .- The blunted cone with flat or 
rounded tip has become increasingly common in hypersonic vehicles . 
Attempts were made to adapt the present method of prediction to the 
general case of truncated cones but with little success. Characteristic 
reflection methods similar to those mentioned by Giese and Bergdolt 
(ref. 29 ) proved to be inadequate primarily because methods that were 
suitable for estimating the overexpansion and subsequent convergence 
toward conical flow at low Mach numbers were completely inadequate at 
high Mach numbers (and vice versa) . Nevertheless one point of possible 
convenience was exposed. So long as the semicone angle of the blunted 
cone does not exceed about 15° -for Mach numbers in the neighborhood of 5 
and below, decreasing to about 8° at Mach numbers near 8, the effect of 
overexpansion in producing a reflex in the shock shape is negligible 
for all practical purposes, and the shock may be calculated by the usual 
procedure with the additional condition that, when its slope deteriorates 
to that which would be produced by the cone in the absence of tip blunt- 
ness, the slope is maintained constant at this value. 
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CONCLUDING REMARKS 


A reexamination has been made of the use of simple concepts for 
predicting the shape and location of detached shock waves . The results 
show that 'simple concepts and modifications of existing methods can 
yield good predictions for many nose shapes and for a wide range of 
Mach numbers . 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., August 28, 1957* 
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APPENDIX A 

DEVELOPMENT OF GENERAL ELLIPTIC EQUATION GIVING THE VARIATION 
IN DETACHMENT DISTANCE BETWEEN 6 Q = 8 det AND S Q = 90° 

FOR CONE-CYLINDERS AND WEDGE-SLABS 


Consider the general elliptic equation 


i (x - l) 2 + 1 (y - k) 2 = 1 
a 2 b 2 


with center at Z,k and semiaxes a and b. 


(Al) 


y 



From the sketch 


k = b + N (A2) 

Differentiating equation (Al) with respect to x and substituting 
equation (A2) into it yields 


dy _ b^ x - l 
dx a 2 y - (b + N) 


(A3) 
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If the slope at x = 0, y = 0 is denoted by q, then frotn equation (A3) 
it follows that 


a 


2 _ 


b 2 l 


q(b + N) 


(A4) 


Also at x = 0, y = 0, substituting equation (A2) into equation (Al) 
yields 


a 2 = 


7 2 b 2 


Equating (A4) and (A5) gives 


b^ - (b + N)' 


b = IqN - IP 
2K - Iq 


(A5) 


(A 6) 


Equation (Al) may be expressed as 




2(b + N) 


y + 


± (x - Z ) 2 - l 1 


b £ 


= 0 (A7) 


which is readily recognized as a quadratic, the solution for which is 
easily obtained. In terms of the shock -detachment parameters, the con- 
stants and variables involved in the solution of equation (A7) are 


x = 5 0 I 5 det 

57-3 


l 

q. 

N 



37-3 S 

-0.5 csc 2 S det 

p( c 90 " l) 


(A8) 
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where p = 0.5 cot 5^ e ^. and C^q is the value of C for 6 q = 90° 
or the flat-face condition; that is, C^q = 0.70 for cone -cylinders 
and = 0.86 for wedge-slabs (two-dimensional). 

Values of x are determined by allowing 5 q to vary between 

and 90° • The corresponding values of y calculated from equation (A7) 
are converted to the detachment distance parameter x’/d 1 by 

fr = p + y (A9) 

d’ 
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APPENDIX B 


METHOD FOR CALCULATING ELLIPTIC VARIATION - 
OF C FOR CUT CYLINDERS AND SPHERES 


Consider first a typical variation of x'/d' for a cone -cylinder 
(or wedge -slat) as 8 q varies between and 90 °, that is, the 

curve BB' in the following sketch: 



As shown previously, the maximum value of x'/d' occurs at point B and 
at this point C = 1. The minimum value of x'/d' occurs air point B' 
and at this point C is equal to that for a flat face or C^q. It 

is readily recognized that C^q and point B 1 also apply to the case of 

the cut sphere (or cylinder) for 0 s 0°, where 9 = 90 ° - 6q, since as 

0 approaches 0° the cut sphere approaches a flat face (limit- of 8 = 0° 
excluded ) } or alternatively, for the case of constant maximum diameter 
D and varying nose radius R, r/d = <» when 0=6° and the face is 
flat . ... 

The value of C for_the complete sphere (that is, C c ^ is always 

less than 1, and the corresponding value of— x'/d' may be denoted on 
the curve BB', for example, at point G. The value of 8 Q corresponding 

to point G is incidental to the development herein, since it merely 
defines the cone-cylinder giving the same detachment - distance as the com- 
plete sphere. The importance of point G lies in the slope at this point. 
For the cut sphere the slope dC/a&Q and the value of C at Sq = 90° 
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are known as is the value of C at 


6^ S 5, 


det ; 


in line with the Busemann 


concept C is assumed to be equal to C c for 8 q ^ S& e t* Thus the only- 

quantity needed to determine the elliptic variation (assumed) in C for 
the cut sphere is the slope dC/d&Q at 5^ e ^_ . No ready solution was 

found for this slope (as contrasted to that for cone-cylinders), but it 
is likely close to that which may be determined from the curve BB' at 
point G. The assumption is therefore made herein that the slope dC/dS^ 


for the cut sphere at 8^^ is equal to that which may be determined 

from this point. The variation of C for the cut sphere in comparison 
with that for the cone -cylinder may thus be sketched: 



c 90 


and the slope at E is equal to the slope at E' . 


The procedure for calculating the elliptic variation is as follows : 
Determine the value of C & at the particular value of M ra from figure 2 

Substitute this value of C c into the following expression and obtain y 


y = p(c c - l) 


(Bl) 


where p = 0.5 cot . Substitute y into the elliptic equation 
(derived in appendix A) and solve for x: 


/V\ v 2 _ 

2(b + N) 

M 4 * 

(b + N) 2 

(b 2 f 

L b2 

y ^ 

b 2 • 


+ -l ( x - ir - 1 

BT 


= 0 


(B2) 
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where, as in appendix A, 

a 2 = - fr 2 * 

q(b + N) 


b _ ZqN - N 2 
2N - 2q 


q = 


-0.5 csc^S 


det 


(B3) 


, , *?.: j d e t 

57.3 

K - p(C 90 - 1) 

and Cgo = 0.86 for the two-dimensional case and C^q = 0.70 for the 
axdsymmetric case. 


Although it is not required in this procedure, if one desires the 
value of 5 q corresponding to the complete - sphere (or cylinder), it 
may he obtained by substituting this calculated value of x into the 
general- relation 


x 


S 0 ~ 5 det 

57-3 


(B4) 


Next, with the value of y from equation (Bl) and the value of x 
from equation (B2), calculate dy/dx from 



(B5) 


This value of dy/dx equals 


which in turn equals 
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And since 


-£C = l iiliz (b6) 

dS 0 P d5 0 

the value of dy/dx from equation (B5) may be substituted into the 
following equation : 


dC = 1 dy 
dB Q p dx 


(B7) 


to obtain dc/d 8 0 which is the desired slope at E and E 1 . (See 
sketch. ) 

It now remains to determine the variation of C with 8 q (or 0 ) . 

In order to do this, the elliptic equation (B2) is used, but note that 
the constants q and N are now 



f (B8) 



and therefore a 2 and b are also different. From equation (b 4) 
and the relation 6 = 90 ° - 9, 


x 


(90° - e) - s flgt 

57-3 


(B9) 


Let 9 vary from 0° to 90° - 6^ e ^., substitute the corresponding values 

of x into the elliptic equation with the indicated change in constants, 
and solve for y. 
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Obtain values of C from 


0 = C c + y (BIO) 

Substituting these values of C into 

fr - c P (Bii) 

gives the variation of x'/d 1 with 8 ^or 5 q) 
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APPENDIX C 


DERIVATION OF METHOD FOR CALCULATING r\ FOR KNOWN SHOCK SHAPES 

The following quantities are known or can "bej determined according 
to the methods suggested in the text: M^, 3> s s , ^det^ an ^" 

Jr = °-5C cot 5 det . 


From reference 9 (see also fig. 9 of this paper) and in terms of 
the present parameters 



(Cl) 


Similarly, the results of reference 9 can be used to obtain 

A 


fa = 


p - 1 (|1 + *9§j) 


p 

3 tan 6 S 


3 ^p 2 tan 2 e s - 1-+ tan q 


(C2) 


Combining equations (Cl) and (C2) yields 



The q uan tities y/d' and x/d' are the ordinates of a point on the 
known shock shape. The shape, however, is usually' given or obtained in 
terms of the maximum diameter D and with distance in the x-direction 
measu red from the ce nter of the spherical Cor circular) nose- or from the 
face - of a flat nose"! This distance will be called P herein^ The 
following simple expressions for obtaining y/d 1 and x/d' from these 
dimensions are given for the spherical (or circular) nose: 
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2- = Z i 

d' D cos & det 


(C4) 


and for the flat"" nose: 


_y 

d* 


y 

D_ 


Values of x/d' 


are obtained from 


x 

d' 



(C5) 


(c6) 


*where the quantity t has the following value for the spherical (or 
circular) nose: 


^ = 


cos 8 


det- 



sin^B 


det 


+ C cos^S, 


det' 


2 s’in • 8 


det 


(07) 


and for the flat nose : 



(C8) 
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APPENDIX D 


PROCEDURE FOR CALCULATING SHOCK SHAPE 


The following quantities are known or can "be determined according 
to the methods suggested in the text: 0, e s , C, q, ^det^ an<a - 

j-j- = 0.5C cot . The value of is obtained from equation (C2) . 

x n 

With the value of — thus determined, it is convenient to determine 
- d f * 

-y- by assigning positive values to the quantity —■ - — from zero 
d d d 

upward. The value of -p- may be calculated from 


y. = i (t d \ 2 

d’ pvU 1 / VdV 


\ 


(Dl) 


In terms of the diameter D the ordinates become 


Hz) 


cos 8, 


det 




(D2) 


and 


* 


, F 
" D 


X 

d* 


x oV.„ « sln26 det + C cos25 det 

^ V *** 2 sin 8 det 


(D3) 


If the slope of the shock at any point is desired, it may be 
obtained from the relation 


dx 


d' 


- ($‘ 


(D4) 


>0 


f wO 


• >v' v 
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Axisymrrietric data 

Serbin, ref. 16 -flat face 
Serbii, ref. 16 -theor.vdue for Mm"® 
Heberle, et at, nef. IO-cone-cylinder 
Heberle, et aL, ref. |0- sphere 
Oliver, ref. 20-sphere 
Cooper and Robinson, Langley I I'HST- 
flat face (unpublished) j 

Moccoll, ref. 35 -cane-cylinder 
llrpubished data, Langley 9"TS - 
flat face 

Sugmoto fevalable in ref. 1 9 -sphere 
Moedtel, ref. 1 1 -flat face, rounded 
comers 

Moectel, ref. 1 1- sphere 

Crawford end McCauley, ref. 38-sphere 



o\ 


4 5 6 

Free-strecm Mach runber, 


Figure 1.- Compilation and general correlation of data on detachment distance for tiro-dimensional 

and axisynmetric nose shapes in air. y = 1.4. 
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Cone or wedge length and detochmetf distance 



Figure 3*- Variation in detachment distance vith semiapex angle for cone-cylinders and vedge- 

slabs at constant Mach number. 



\J4 

oo 


NACA TN 4170 


Detachment distance, 


NACA TN 4170 


39 



Free-stream Mach number, Mjq 


Figure 4.- Predicted variation in detachment distance with Mach number 
for cone- cylinders at constant semiapex angle. 
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Prediction of Kim's results (ref. 24) for cut circular cylinder (two-dimensional) at 

H, = k. 
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Figure 8.- Prediction of results at = 3*35 for a sphere that is 

effectively cut. (Diameter of actual models held constant and radius 
of nose: varied. ) 




Figure 9 • - General features of the method of Moeckel (ref. 9 ) ado pted in present analysis 
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(a) Sphere. 


(b) Circular cylinder (two-dimensional) . 


Figure 10.- Comparison of variation of tj with Mach number obtained in present analysis with 

those obtained in other analyses. 
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Figure 11.- Effect. of C upon shock shape and the possible ambiguity 
that can occur with the present— analysis when the 6hock is presented 
in this form. ' M™ = 1.94. ; 
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(a) Hemispherical nose. 

16 " R ' ediCtl0n ° f ^ ^ - looat lon for exiBymmetric noses at - . 6.8. 
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Figure 17.- Prediction of shock 6 hape and location for hemisphere-cylinder at = 7.7. 
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